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Abstract 

Spatial localization of the electrons of an atom or molecule is studied in models of non- 



relativistic matter coupled to quantized radiation. We give two definitions of the ionization 

>• ■ 

threshold. One in terms of spectral data of cluster Hamiltonians, and one in terms of 
CN . minimal energies of non-localized states. We show that these two definitions agree, and that 

' the electrons described by a state with energy below the ionization threshold are localized 

in a small neighborhood of the nuclei with a probability that approaches 1 exponentially 
fast with increasing radius of the neighborhood. The latter result is derived from a new, 
general result on exponential decay tailored to fit our problem, but applicable to many 
non-relativistic quantum systems outside quantum electrodynamics as well. 

a 

1 Introduction 

If an atom or molecule is in a state with total energy below the ionization threshold, then all 
electrons are well localized near the nuclei. In non-relativistic quantum mechanics this finds 
its mathematical expression in the discreteness of the energy spectrum below the ionization 
threshold and in the exponential decay of the corresponding eigenfunctions. When the electrons 
are coupled to the quantized radiation field, then there is no discrete spectrum anymore and the 
ground state is the only stationary state [HE]. Nevertheless, all states in the spectral subspace 
of energies below the ionization threshold are exponentially well localized as functions of the 
electron coordinates. To prove this is the main purpose of this paper. Localization of the 
electrons below the ionization threshold is necessary to justify the dipole approximation [2], 
and it plays an important role in proving existence of a ground state [21 E3 E| an d Ior Rayleigh 
scattering [5]. 

The ionization threshold is the least energy that an atom or molecule can achieve in a state 
where one or more electrons have been moved "infinitely far away" from the nuclei. To give 
*Work partially supported by U.S. National Science Foundation grant DMS 01-00160. 
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a more precise definition we need a mathematical model for atoms and molecules. A (pure) 
state of N electrons and an arbitrary number of transversal photons shall be described by a 
vector in the Hilbert space TLn = TLei ® where TL e i is the antisymmetric tensor product of 
N copies of L 2 (R 3 ;C 2 ), appropriate for N spin-1/2 fermions, and J- is the bosonic Fock space 
over L 2 (M 3 , C 2 ; dk). The nuclei are static, point-like particles without spin. Let Hn denote the 
Hamilton operator generating the time evolution in TLn, and let H N be the same Hamiltonian 
without external potentials (nuclei). We assume that the dynamics of the electrons is non- 
relativistic and that the forces between material particles (electrons and nuclei) drop off to 
zero with increasing distance. In view of the latter assumption a natural definition for the 
ionization threshold t(Hn) is 

t{H n ) : = min {Ejy—N' + E^>}, (1) 

where E^-n' = inf c(-Hjv'-jv')j E n , = m ^ a (-^N')' anc ^ ^N=o = 0. Let m be the mass of the 
electron and let \x\ = (X)j=i ^f) 1 ^ 2 f° r x £ We prove that, for all real numbers A and f3 
with A + (3 2 /{2m) < t(H n ), 



< oo, (2) 



and that in states with energy above t{Hn) the electrons will not be localized in general. 
Thus t{Hn) is in fact a threshold energy separating localized from non-localized states. The 
question of whether the binding energy t{Hn) — En is positive or not, is not addressed in this 
paper, see however [7]- 

Our proof of (J2J) consists of two independent parts. First we give an alternative definition 
of the ionization threshold which better captures the idea of a localization threshold, and we 
prove exponential decay below it. Then we show that the two definitions agree. 

The alternative definition is as follows. Let Dr = {ip £ D(H)\ <p(x) = 0, if \x\ < R}, and 
define a threshold energy T,(Hn) by 



E(£fr) = lim inf (<p,H NV >) . (3) 

R— +oo \ipeD R , ||v||=l / 

Derealization above S(iJjv) is obvious, and localization below T,(Hn) will be derived from 
the only assumptions that Hn is self-adjoint, bounded from below, and that 

[[H N ,f],f] = -2\Vf\ 2 (4) 

for all bounded smooth functions f(x) with bounded first derivatives. The latter assumption is 
satisfied for the positive Laplacian (—A), and hence for all operators — A + I with [[/, /], /] = 0. 
Examples include the commonly traded models of non-relativistic atoms coupled to quantized 
radiation, as well as many Schrodinger operators outside quantum electrodynamics. 
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The second part of the proof, that t(Hn) = T,(Hn), is the hard part. The inequality 
t{Hn) < £(-ffjv) requires localizing both the electrons and the photons, and in particular 
their field energy. This was done in FF]. To show that t(Hn) > S(-f/jv) we construct suitable 
(compactly supported) minimizers (fo and ip^ of Hn-n' and H N ,, respectively, where (p^ 
is localized at a distance R from the origin. We than merge these states into a single state 
tpn S TLn- The problem is to do this in such a way that (^r, Hn^r) = (pq, Hn-n /( Po) + 
+ o(l) asR^oo. 

In the context of QED the first result of the form © is due to Bach, Frohlich and Sigal 0, 
who proved exponential binding for small coupling and away from the ionization threshold of 
Hn with zero coupling. The threshold energy t{Hn) was introduced in [7j where it was shown 
that En is an eigenvalue of Hn if t(Hn) > En- The paper [7] also contains an easy argument 
showing that eigenvectors of Hn with eigenvalues below t(Hn) exhibit the exponential decay 
implied by (J2J. For iV-particle Schrodinger operators the ionization threshold defined by the 
analog of (^TJ is the least point of the essential spectrum. This is known as the HVZ-Theorem 
[S]. That the analog of © also characterizes the beginning of the essential spectrum is a result 
due to Arne Persson ^U]. Exponential decay for TV-body eigenfunctions with discrete energy 
was first proved by O'Conner j^]. See Agmon's book 1 for more results on the exponential 
decay of solutions of second order elliptic equations. 

Section [2 contains the general theorem on exponential decay in an abstract Hilbert space 
setting. In Section |21 this result is applied to quantum electrodynamics and the main result on 
equality of the thresholds is formulated. Its proof is given in Section The Appendix collects 
technical results and notations used in the proofs. 



2 The Abstract Argument 

In this section q : D x D — > C denotes a densely defined, closable, quadratic form that is 
bounded from below and defined on a domain D C H in a Hilbert space H. We assume that 
TC is a closed subspace of a Hilbert space L 2 (M n ) ® T and that TL is invariant with respect to 
multiplication with bounded (measurable) functions that depend on x € W 1 , only. Here 
J- is an arbitrary, additional Hilbert space. In our applications T will be the tensor product 
of spin and Fock space and 7i the subspace with the symmetry required by the nature of the 
particles. 

On the quadratic form q we make the further assumption, that for each / £ C°°(M n ;M) 
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with /, V/ G L°°(M n ) and with f(x) = there exist constants a and b such that 

(i) /DCD 

(**) /¥>)l < ag(^, + & 

(iii) q(f<p, if) + gfo /V) - 2g(/y>, jV) = -2 (p, | V/| V> 

for all (/? G D. Requirements (i) and (ii) are mild technical assumptions which ensure that 
property (iii) extends to all ip in the domain of the closure of q. Equation (hi) is the basis of 
the so called IMS (localization) formula for Schrodinger operators. To verify it for a quadratic 
form q that is defined by a symmetric operator H : D C TL — > TL it is useful to know that 
f 2 H + Hf 2 - 2fHf = [[H, /],/]. Assumption (iii) then becomes 

[[HJ],f] = -2\Vf\ 2 , 

which holds for the positive Laplacian (—A) and hence for all operators —A + I in TL with 
[[/, /], /] = 0. Some examples, other than those in the next section, are H = (— iV + A(x)) 2 + 
V(x) with a classical vector potential A(x) and scalar potential V(x) (choose T = C), and 
Schrodinger operators with restricted domains £1 C 1" (TL = L 2 (Q) C L 2 (W n ) C), or with 
potentials that are constant away from a strip, as in wave guides defined by potential wells. 
Given R > let Dr = {ip G D : ip(x) = for \x\ < R} and define 

= inf q(cp,<p) and £= lim (5) 
<p£Dn, ||^||=1 R^oo 

The numbers T,r are finite because q is bounded from below and because, by (i), Dr is not 

empty. But E may take on the value +oo. 

Theorem 1 (Exponential decay). Suppose the quadratic form q introduced above satisfies 
the assumptions (i), (ii), and (Hi), and let H denote the unique self-adjoint operator associated 
with the closure of the form q. If X and (5 are real numbers with A + [3 2 < S, then 

3^E X (H) 

where E\(H) is the resolution of the identity for H . 

Remarks. (1) For Schrodinger operators —A + V on open domains C W 1 with Dirichlet 
boundary conditions and with V- <C —A the above theorem implies that the spectrum below 
S is discrete. In fact (—A + l)~ 1 / 2 e~^l x l is compact and hence so is E\(H) = E\(H)(—A. + 
l)i/2 (_ A + -q-i/2 e -/3M e P\x\E x (H) for A + /3 2 < S. 

(2) Everything in this section holds equally for any norm |x| on W 1 that is induced by an 
inner product x ■ y, provided that A is used to denote the Laplace-Beltrami operator with 
respect to the metric g(x, y) = x ■ y. 

The following proof is inspired by the proof of binding in Bach et al.|2]. 



< oo, 
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Proof. Let Q(H) C TL denote the form domain of H, i.e., the domain of the closure of q. We 
use q to denote the closure of q as well. Q(H) is the closure of D with respect to the form norm 
|| • || q associated with q. By assumptions (i) and (ii), multiplication with a bounded function 
/ G C°°(R n ) with bounded derivatives is a bounded linear operator on (D, \\ ■ \\ q ) and hence 
extends to a bounded linear operator on (Q(H), \\ ■ \\ q ). In particular 

fQ(H) c Q(H) (6) 

and (iii) extends from D to Q(H). 

Let E = 'mfcr(H). We may assume E > E, for otherwise there is nothing to prove. Let 
X2R denote the characteristic function of the set {x G W 1 : \x\ < 2R}. We first show that 

H R :=H+ (E h - E) X 2R >^R~Jp ( 7 ) 

for all R with T,r > E and some constant C. Pick j 1 ,j 2 S C°°(R + ) with j\ + jf = 1, 
supp(ji) C {i < 2} and supp(j2) C {t > 1}. Let = jj(|x|/-R). Then by © and since 

(iii) holds on Q(H), 



i=l 

2 2 

= ^3i,RHR3i,R - ^2 \ v kn 



(8) 



1 = 1 1=1 

in the sense of forms on Q(H). By definition of E# and by the construction of j^r, 

h,R H R ji,R = 3i,r( h + ^R~ E )ji,R > ^RJl,R 
32,rHr32,R > 32,RHj2,R > ^RJ2,R- 

Hence (JZJ) follows from (jSJ) and from (Vj'^i?! = 0{R^ 1 ). 

Let A := [inf a(H), A], where A+/3 2 < E, and pick R G R so large that A+/3 2 < Y, R -C/R 2 . 
This i? is kept fixed in the following. Let S := T, R — C j R 2 — (3 2 — A > 0, and choose a function 
SA G C£°(R;[0, 1]) such that g A = 1 on A and supp(#A) C (-oo, A + 5/2]. Then, by Q, 
9a{Hr) = and therefore 

g A (H) = g A (H) - g A (H R ). (9) 
We now show that e^' ($a(-£0 — g A (H R )) is bounded. To this end, we define 

and show that, e* (g A (H) — g A (H R )) is bounded uniformly in e > 0. Note that / G C°°(R ra ), 
is bounded and that | V/| < f3. Let g A be the almost analytic extension g A (x + iy) = (g A (x) + 
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iyg' A (x))'y(y) where 7 £ Cq°(R) equals one in a neighborhood of y = 0. By the almost analytic 
functional calculus (see 03) 

g A (H)=-~ J ^(z-HyUxdy 
and hence, using Q and a resolvent identity we can write 

e f g A {H) = ^ J ^e f (z-H R )- 1 e- f e f (^ R -E) X 2 R (z-H)- 1 dxdy 
whose norm we estimate from above as 

\\efg A (H)\\ < sup \\ef(z-H R y 1 e~f\\\\ef X 2R\\ 00 (Z R -E) 

zesupp(g) 

(z — H)" 1 ]] dxdy. 



x — 

7T 





dg 


/ 


dz 



The norm \\e* X2i?||oo is bounded uniformly in e > and the integral is finite. To estimate 
He* (z — H R )~ 1 e~f\\ let H R j := e* H R e~f with domain D(H R j) = e* D(H) and note that 



(z-Hrj)- 1 =ef(z-H R )- 1 e~f 

as can easily be seen by direct computation. In particular, the resolvent sets p(H R j) and 
p{H R ) coincide. Let (p £ D(H R j) C Q(H) and ||</>|| = 1. Then 



2Re (p,H RJ p) = (p,(e f H R e f + e f H R e f )p 

= ^ip,e~ f (e 2f H R + H R e 2f )e~ f ip 
= 2(p,(H R -\Vf\ 2 )p) 

where (in) was used in the last equation. In conjunction with (J2J) this shows that, for z S 
supp(<7), 

Re {<p, {H RJ - z)tp) >Z R - C/R 2 - (5 2 - Re(z) > 5/2 (10) 
and hence that \\(H R j — z)tp\\ > <5/2||(£>||. Since p(H R j) = p{H R ) D supp(g), it follows that 

\\(z-H RJ )- l \\<2/5 

for z S supp(^), which completes the proof. □ 

3 Atoms Coupled to Quantized Radiation 

In this section we apply the abstract result of the previous section to systems of N charged, 
non-relativistic quantum particles, interacting with the quantized radiation field. Since we are 
mainly interested in the case of electrons in the field of static nuclei, the bulk of the exposition 
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deals with this case. At the end we comment on the more general case of particles from different 
species. 

In the "standard model" of non-relativistic QED the Hilbert space of a system of N elec- 
trons and an arbitrary number of transversal photons is the tensor product 

H N = A^L 1 L 2 (M 3 ; C 2 ) <8> Tj 

of the antisymmetric product of N copies of L 2 (R 3 ; C 2 ), appropriate for N spin-1/2 fermions, 
and the bosonic Fock space Tf = © n >o <8>" L 2 (R 3 , dk; C 2 ), where the factor C 2 accounts for 
the two possible polarizations of the transversal photons. Let T>n C Hn be the subspace of 
sequences ip = (<po, <p±, . . .) where 

y n € C^ a ((R 3 x {1,2}) N ;C) CL 2 (M 3 ,C 2 ) 

and ip n = for all but finitely many n. The index a indicates that the functions are anti- 
symmetric with respect to permutations of the TV" arguments, and Lq(R 3 ;C 2 ) is the space of 
compactly supported L 2 -functions. Clearly T>n is a dense subspace of Hn- 
The Hamilton operator Hn : T>n C Hn —> Hn of our system is given by 

N 

Hn = J2(Pj + V^A{ Xj )f + | V5o-j • B{ X j) + V + H f , (11) 

3=1 

where pj = —iV Xj , A(xj) is the quantized vector potential in Coulomb gauge evaluated at the 
point Xj, B(xj) = curlA(xj) is the magnetic field, <jj the triple of Pauli matrices (a^~\a^\ <?^) 
acting on the spin degrees of freedom of the jth particle, V is a real-valued potential, and Hf 
is the Hamilton operator of the field energy. The parameter a is the fine structure constant 
and the coupling constant g 6 R is arbitrary, to allow for a simultaneous treatment of the 
interesting cases g = 2 and g = 0. 
Formally A(x) is given by 



A(x) = W -^=e x (k) \e ik -*a x (k) + e-**aj(*) 

> , n J\k\<A y/\k\ 



A=l,2 



d 3 k 



where A < oo is an arbitrary but fixed ultraviolett cutoff. For every k ^ the two polarization 
vectors £\(k) £ R 3 , A = 1, 2 are normalized, orthogonal to k and to each other. 

The operators a\(k) and at (k) are the usual annihilation and creation operators, satisfying 
the canonical commutation relations 

[a x (k),a;(k')\ = 5^5(k - k'), [a*{k),a*{k')\ = 0. 
In terms of a\(k) and a* x (k) the field Hamiltonian is given by 

H f =J2f d 3 k\k\al(k)a x (k). 

X=l,2 
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See Appendix [B] for mathematically more proper definitions of A(x) and Hf. 

The potential V is the sum of the external potential and the Coulomb two-body potentials 
for each pair of electrons. For the purpose of the results to be proved in this section, however, 
it suffices to assume that 

(HI) V G Lf oc (R 3N ; R), and V_ < e(-A) + C £ for all e > 0, 

and, of course, that V is symmetric with respect to permutations of the particle coordinates. 
The Hamiltonian Hn is a symmetric, densely defined operator and by Lemma|HJ it is bounded 
from below. The quadratic form q(<p, ip) = (ip, HjfiftS with domain D = Vjq is therefore 
bounded below and closable and hence the theory of the previous section applies, once we have 
verified assumptions (i), (ii), and (iii). The unique self-adjoint operator Hn associated with 
the closure of the quadratic form q is the Friedrichs' extension of fljy. The thresholds and 
S associated with the form q are now given by 



Er(H n ) = inf (ip,H N tp 

<peT>N,R, IMI=i N 

E(H N ) = lim E R (H N ) 

R—rOO 

where 2?7v 5 _r := {<p 6 Pat : (p(X) = if \X\ < R}. The following theorem is a corollary of 
Theorem ^ 

Theorem 2 (Exponential decay in QED). Assume Hypothesis (HI) is satisfied and let 
Hn be the Friedrichs' extension of the symmetric operator Hn : P/v C T~Ln 7~Cn given by 
Eq. If X and (5 are real numbers with A + /5 2 < J^(H]y), then 



e^ x \E x (H N ) 



< oo. 



Proof. It suffices to verify the assumptions (i), (ii) and (iii) in the previous section. Suppose 
/ G C°°(lR 3Ar ) with /, V/ G L°°(R 3Ar ), and f(X) = f(\X\). Then fV N C V N is obvious from 
the definition of T>n- Property (ii) follows from 



f{pi + V^A{x t )ff < 2\\f\\l (p i + V^A(x t )) 2 + 2\\V x J\\ 2 00 
fH f f < Wff^Hf 
fVf < H/foo^t-, 

from Lemma [7| and Lemma |S1 The proof of (iii) is a straightforward computation using that 
f 2 H N + H N f 2 -2fH N f=[[H N J],f}. □ 

Remark. The above theorem and its proof can easily be generalized to systems of N 
particles from n < N species, with different masses m;, charges, and spins. Theorem |21 then 
equally holds with the new norm \X\ = (^2iLi 2mixf) 1 ^ 2 in the factor e^ x V 
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Our next goal is to establish a relation between S(i?jv) an d spectral data of cluster Hamil- 
tonians. To this end, we impose the following additional assumption on V: 



(H2) 



v ( x ) = Eili v ( x i) + Ei<j w ( x i ~ x j) where v,w £ L 2 oc 
and lim\ x \^ 00 v(x) = 0, lim^^wix) = 0. 



If the external potential v is associated with a particle sitting at the origin x = 0, then 
these assumptions can be understood as saying that the interaction between spatially separated 
clusters of particles drops off to zero as the inter-cluster distance increases to infinity. The 
limitation to two-body forces in (H2) is not necessary. 

Theorem 3 (Equivalence of ionization thresholds). Assume (HI) and (H2), and let 
En_ni = inf ct(Hn_n/), and E^ = 'mia(H^) where all external potentials are dropped in H^. 
Then 

E(H N ) = mm{E N - N , + E° N ,}. 

N'>1 

The proof requires, in particular, localizing the field energy in neighborhoods of the elec- 
trons. In order to control the localization errors which thereby arise we need an infrared cutoff 
in the interaction. That is, we first prove the above theorem in the case where all interactions 
of electrons with photons of energy less than an arbitrary small, but positive constant fi have 
been dropped from H^. The theorem then follows in the limit fi — ► 0. 



4 IR-Cutoff Hamiltonians 

In this section we prove Theorem |31 by first establishing an analogous results for Hamiltonians 
with an infrared cutoff [i in the interaction. Theorem then follows in the limit fi — ► 0. 

The infrared cutoff Hamiltonians Hn a, H > 0, are defined in the same way as with the 
only difference that the vector potential A(x) and the magnetic field B(x) in Hn are replaced 
by 

^ r 1 r ., 

d 6 k 



M*)= >J / 4^ £ a(£0 e^ x a x (k) + e-^ x a* x (k) 

//i<|fc|<A 



E 

A=l,2 



and Bn{x) = curL4 M (x). To separate the soft, non-interacting photons from the interacting 
ones we use that Tf is isomorphic to T{ ® T s where Ti and T s denote the bosonic Fock spaces 
over L 2 (\k\ > fj) and L 2 (\k\ < fi) respectively. Let Hi = A^L^IR 3 ; C 2 ) Then the 

Hamilton operator can be written as 

H N ^ = Hi®l + l®H s f on H = H i ®J : s (12) 

if we identify T with Ti® T s . Let T s n denote the n-boson subspace of T s and let Q s be the 
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vacuum of T s . Then l|12[l and the positivity of Hj = Hf\T s imply that 

infer {Hn,ii) = hif (mia(HN,n\'Hi® F s ,n)) 

n>0 ' (13) 

= mfa(H N jHi®[n s ]), 

where [Q s ] is the space spanned by £l s . This will allow us to drop the soft bosons in all 
approximate energy minimizers. 

Lemma 4. There exists a constant CV,a> depending on A,N,g and a, such that 
±(H N -H N ^)<fi 1 / 2 C NA i^p? + #/ + lJ for 0<fi<l. 
Proof. By definition of Hjy and Hjv>, 

N 

H N -H N}tl = Yj 2 Vapi ■ {A( Xi ) - A^x,)) + a(A( Xi ) - A^(xi)) ■ (A( Xi ) + A^Xi)) 



i=i 
9 
2 



where we used that A(x) and ^(x) commute. The differences A{x) — A^(x) and B(x) — B^(x) 
can be seen as a vector potential and a magnetic field with an ultraviolett cutoff fj,. Hence the 
lemma follows from Lemma with A = /i and e = /U 1 / 2 . □ 

Lemma 5. (%) < oo if and only ifT,(H^) < oo and in i/iis case t/iere exists a constant 

Ca depending on the parameters A, N, a, and g, such that 

|E(fliv >M ) - £(#jv)| < Cam 1/2 , »/a* < 1- 
(ii) There exists a constant Ca depending on the parameters A,N,a, and g, such that 

\t(H n ^) - t(H n )\ < C A ^ 2 , i/m < 1. 
Proof. By Lemma0]and Lemma[SJ there exist constants C and D, independent of fi, such that 

H n ^<H n + ij}I 2 {CH n + D) for /x < 1. 

It follows that 

V{H Njlx ) < E(H N ) + ^l 2 {C^{H N ) +D) for fi < 1 

and, in particular, that E(iJjv,/Lt) < oo if T,(Hn) < oo. Since the roles of i?jv> and i5/jv are 
interchangeable, (i) follows. The proof of (ii) is similar. □ 

Theorem 6. Suppose assumptions (HI) and (H2) on V are satisfied. Then 

T,(H N ^) = t(H n ^) for all fi > 0. 
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In conjunction with Lemma [51 this theorem proves Theorem |31 

Proof of E(iljv"u) > t {Hn4l)- The key element for this proof is Theorem^ whose long proof is 
given in ,7]. Here we merely show how Y,(Hn^) > t(Hn^) follows from Theorem HO We may 
certainly assume that T,(Hjy^) < oo. By the argument ()13j) we may restrict to Tii <S> [Q s ] 
for the computation of ^^(Hn^). By Lemma|H] 

N f < -H f < -(2H N „ + D) on Hi ® [Q a ] 
and hence by Theorem |SJ 

H N ^>r(H N ^)-o(R°)(H N ^ + C) on v N>R n(Hi®[n s ]). 

It follows that 

^r(H n ^) > r(H Ntlt ) - o{R°)(E R (H N ^) + C) 
and the desired result is obtained in the limit R — > oo. □ 

An important role in the following proof is played by the identification operator / : T®T — * 
T which collects all photons in the first and second factor of J 7 ® !F, and gathers them in a 
single Fock space. For the precise definition of I, and for notations in the following proof that 
have not yet been introduced, see Appendix iBl 

Proof o/ S(i7jv )At ) < t(Hn^). In the following the subindex u is dropped. We need to show 
that 

lim Tir(Hn) < -Ejv-Tv + E N , 

R^oo 

for all N' > 1. The strategy is as follows. First we construct approximate minimizers ipo and 
ifoo of Hn_.ni and H N , respectively, with the property that the electrons and the photons 
described by (po and ifoo are compactly supported. Then, by a translation (p^ — > Tjnp^ of 
both the electrons and the photons in 9300 we may achieve (ignoring the Pauli principle) that 

ip R = I(<p <8> Tr^oo) G T>n,r, and \\ip R \\ = 1, 

where Tr^qo is still an approximate minimizer of H N , by the translation invariance of this 
Hamiltonian. 

Second we show that 

(ip R , H N ip R ) < (ipo, Fjv-at'Vo) + (<Poo, Hn'Voc) + o(R°) R -> 00 

which concludes the proof. To incorporate the Pauli principle one needs to anti-symmetrize 
I(ipo (8> T R (poo) with respect to the N electron variables (xj, Sj) G K 3 x {1, 2}, i = 1, . . . , TV. 
After normalization, this will lead to the same value for the energy {ip R , H^^Pr) as without 
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anti-symmetrization, because the electrons in ipo and T^ip^ are disjointly supported and the 
Hamiltonian is local. Therefore we don't need to anti-symmetrize. 

Let e > be given and fixed in the following three steps, and let y denote the position 
operator y = iS/^ in the one-photon Hilbert space. For simplicity the irrelevant parameters a 
and g are dropped henceforth. 

Step 1. Given e > there are normalized states ipo £ Vn-n 1 and <^oo S V^i such that 

(i) (^o, H N _ N ,ip ) < E N _ N , + e/2 and (^oo, #j\r/<Ax>) < e n' + e / 2 - 

(ii) Both (ipo, Nfipo) and (</3oo, -W/Voo) are finite and bounded by a constant that is indepen- 
dent of e > 0. 

(iii) (/?o and i^oo have compact support as functions of the electronic configurations Xn-n> 6 
R 3(N~N') and XAr , G R 3iV' reS p ec ti V ely. 

(iv) There exists an Rq such that 

where XRa 15 the characteristic function of the ball {y 6 M 3 : < i?o}- 

Proof of Step 1. The properties of the Hamiltonians that are relevant, are shared by -H^v-tv 
and Hpji- So it suffices to prove existence of ipo . Let Ho := H^-n 1 and Eq := En-n' fo r short. 
Let xp be the operator of multiplication with x(\X\/P) on Hn-n' where % £ C°°(^+)) 
x(i) = 1 for t < 1, = for i > 2 and < x < 1. Let jr be the operator of multiplication 
with x(\y\/R) on L 2 (M. S , dk). Existence of </?o with property (i) and (ii) follows from the fact 
that Pjv^jy/ is a form core of Hq, the argument and Lemma El If we now show that 

(xp<Po, (#o - E )xp(po) — ► (ip ,(H - Eo)ifo) (14) 
(r(ip)XP¥ :, o, (#o - ^o)r(jR)xPVo) {xp'Po, ( H o - E )xp 1 Po} (15) 

then (iii), and (iv) will follow, because, by the strong convergence xp ~* 1 and T(Jr) —* 1 the 
norm ||r(jp)xp<A)|| is close to 1 for large P and large R. 
Properties (|14[) and l|15|l follow from 

lim [H ,xp]vo = (16) 

P^oo 

lim (iV / + l)- 1 /2[F ,r(i fl )] X P^o = (17) 

(to be proven shortly) by commuting the operators xp and r(jpj through Hq—Eq and using (ii) 
and that s — limp^oo xp = 1 and s — limp^oo T(jp) 2 = 1. Note that T^RjXpT^N-N' C V^-N'- 
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Equation (|16|) follows from 

N-N' 

[Ho, Xp] = Yj (-^)V^Xp • (Pj + A{xj)) - A XjX P 
i=i 

using V xj xp = 0(P _1 ), A Xj xp = 0(P -2 ) and Lemma|Hl 
To prove (|17|1 we write the commutator as 

N—N' 

[Ho,T(Jr)] = Y, { (ft + A{xi))[A(xi),T(j R )\ + [A( Xi ),T(j R )](pi + A{ Xi )) 

1^1 1 (18) 

+ [a i -B(x i ),T(j R )] + [H f ,T(j R )]] 

Using that H f = dT(\k\), the last term in (JTHJ) restricted to (g>™L 2 (M 3 ) is given by [H f ,T(j R )} = 
Ya=i • -®[l^l> 3r\ ■• -®jRi the commutator being the Zth factor. Since ||[|A;|, = 0(i? _1 ) 
it follows that \\(N + l)- 1 /2[ i 7 /?r ( J - jR )]( A r + 1)" 1/2 || = 0(P _1 ), and hence, by (ii), that the 
contribution due to Hf is of order R . To deal with the first two terms in (JT8|) note that, by 
(J2U) and (EH), 

MarO.rfo)] = a*((l-j R )G Xi )T(j R ) - T(j R )a((l - j R )G Xi ) 

where 

lk fl (C/fl - l)«**)XJ'(iV r + 1)~ 1/2 || < sup \\(j R - l)G Xi \\ ^0, as P^oc. 

|xi|<2P 

It follows that the terms in (|18|) which are quadratic in give vanishing contributions, as 

the factors A(x{) outside the commutators can be controlled by (N + l) -1 / 2 . To show that the 
terms in ()18|) with an operator pi vanish in the limit R — » oo it suffices to add to the above 
arguments that pi[A(xi),T(j R )] = [A(xi),T(j R )]pi because pi commutes with A{xi) and T(j R ), 
that PiXP = XPPi ~ i^iXp an d that ||j>i</?o|| < oo by Lemma |H1 The term involving B(xi) is 
dealt with similarly. 

Step 2. Let e, ipo, and be as in Step 1. Pick Rq so large that, with xr is as in Step 
1 (iv), tp = r(xn )vo, Woo = r(xi? )tAx>, <Po(Xn-n>) = if \X N _ N >\ > Rq and tpoo(X N i) = 
if \Xn'\ > Ro- Let R> Ro and pick a vector d £ R 3 with |d| =3. Let Tr : TLn-n 1 — ► 7~Ln-N' 
be the translation 

Ti? = exp f-sPd • { + P/}j 

where P/ = dr(fc) is the total momentum operator of the photons. Then 

ip R := I((p (g> T R <poo) £ V N)2R . 
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Step 3. If R > R then \\i/> R \\ = 1 and 

(tpR, H N ip R ) = (<p , H N ^ N np ) + ((foo, H%,Lpoo) + o(R°), R -> oo. 

In particular T,r(Hn) < E^-N' + E 1 ^, + 2e for all R, which proves the theorem. 

Proof of Step 3. By construction of (po and TRip^ the photons in these states have disjoint 
support if R > Ro . Therefore 

(^r,^r) = (I(<Po 8)T/ ? ^ oo ),/(v3 ® Tr^oo)) 

= (<p TrIPoo, ip TRipoo) 
= (<P0, V0> (<foo, foo) = 1- 

In the following this property of /, that it acts like an isometry on product states with photons 
supported in {|y| < Rq} and {\y — Rd\ < Rq} respectively, will be used repeatedly and tacitly. 
Writing H f = £ A=1 2 / \k\a* x (k)a x (k)d 3 k and using (|2*3)) one gets 

(ipR,H f Tp R } = ((p ,Hfipo) + {TnpocHfTjitpoo) 

+2 Re / \k\ {a x (k)ipo,(p ) {ifoo^xik)^) e imk d 3 k 

A=l,2 

where T^a(k)TR = e lRd ' k a(k) was also used. The third term converges to zero as R — > oo by 
the Riemann-Lebesgue lemma, because the integrand is in L 1 (M 3 ,C 2 ). 

Since the distance of the electrons described by Tr^p^ to the origin and to the electrons in 
tps is bounded below by 3R — 3Ro, we have, by assumption (H2), that 

(4>R, V N 1p R ) = ((p , V N -N ,( Po) + ^ (^Voo, W ( x i ~ x 3) T RVoo) + o(R°), (R -> oo), 

i<j 

as desired. Next we compare 

N 

Yj (Pi + A (xj)?^R) 

3=1 

with 

Y {^,{Pj + A{x j )) 2 ip G ) + Y ( T -R^oo, (Pj + A{x j )) 2 T R p OQ ) . 

j<N-N' j>N-N' 

To this end we write A(xj) = a(G x ) + a*{G Xj ) and use that 

fa + A( Xj )) 2 = pj + 2 Pj ■ a(G Xj ) + 2a*{G Xj ) ■ Pj 



c .) 2 + a*{G Xj f + 2a*(G Xj )a(G Xj ) + \\G Xj \\ 2 . 
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Let j < N — N', then using and again disjointness of the supports of the photons in ipo 
and Tr^qo, one finds that 

{i>R, {Pj + A {xj)f^R) = (tpo, {Pj + A(Xj)fip ) 

+2(<p ®T R (p oo ,pjtp ®a(G X:j )T R (p oo ) + h.c. 
+2(<p ®T R (p 00 ,a(G X: .)cpo® aiGx^Tmpoo) + h.c. 
+ {<po®TR(p 00 ,p ®a(G x .) 2 T R (p 00 ) + h.c. 
+ (a(G x .)(p ®T R <p oo ,(pa®a(G x .)T R (p oo ) + h.c. 
+ (v? ® alGx^TRtpoo, ipo tg> afix^TRtpoo) . 

All terms except the first one vanish in the limit i? — > oo. In fact, 

aiGx^TRPoo = T R a(G :ri _ M )r(x/? )v?oo 

= T R r(xRo) a {XRo G Xj-Rd)<Poc, 

and since < i?o if yoC^ij • • • i x N-N') 7^ 0, we can multiply this in all the above terms with 
XRo( x j)- But then, by ((20)) and using the notation G\{k) = \k\~ l I 2 E\{k)xk{k) 

\\XR (xMXR G Xl -Rd)(N f + 1)- 1 / 2 !! 2 

< sup \G x (xj - Rd-y)\ 2 dy -> (R oo). (19) 

The case where j > N — N' is dealt with similarly. The only difference there is that 
\xj — Rd\ < Ro in the support of TRip^ and the photons in tpo have support in \y\ < Rq. Hence 
(|19|) will be replaced by 

WXRoixj - Rd)a(xR Gx j )(N f + l)-^ 2 \\ 2 

< sup Yl [ \G x { Xj -y)\ 2 dy^0 (R -» oo). 

The terms involving B{xi) are dealt with similarly. □ 

A Important Estimates 

Lemma 7. For all A > 0, e > and a// i£l 3 , 

A(x) 2 < 327rA(iT/ + A/4), 
ia-S(x) < eH f + ^-A 3 . 

For the proof see 7 . This lemma holds equally for A^(x) and B fJi (x) with > 0. 
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Lemma 8. Let C = 1 + 32vraiVA and D = 8iraNA. Then, for all y, > 0, 

TV r N ^| 

J>? < C \ Y,{Pi + V^A^Xi)) 2 + H f \+D. 
i=l U=l J 

Furthermore, if V- < ep 2 + C e for all e > 0, then there exist constants D(e), depending on 
a, g, N, A and e, but not on y, such that 

N - 1 

£(pi + V^A^Xi)) 2 + V+ + H f \ < (1 + + 15(e). 

i=l J 

Proof. The first part follows from < 2(pi + y / aA M (xj)) 2 + laA^ixi) 2 and Lemma The 
second bound follows from the first and Lemma D 

Theorem 9. Suppose the negative parts V- and w- of the external potential v and the two- 
particle interaction w as functions in M 3 drop off to zero as \x\ — > oo. Then for all values of 
the parameters N, A, a, g and \i > 0, there exists a functions f(R) and a constant C, depending 
on these parameters, such that 

Hn^ > r(H Ntli ) - f(R)(H Ntli + N f + C) on V NjR 

where limn-^ f(R) = 0. Here t{H n ^) = inf/v'>i [inf cr(H N _ N < :fl ) + inf a(H%^)]. 

This theorem is a variant of Corollary A. 2 in [7], where we used the positivity of the photon 
mass to estimate Nf in terms of H^. Thus the error term in Corollary A. 2 of [7] depends on 
the photon mass. This was overlooked in |7J leaving a gap in the proof. Theorem |3] combined 
with Theorem 5.1 in [7j closes the gap. 



B Fock Space and Second Quantization 

Let f) be a complex Hilbert space, and let (g)™f) denote the symmetric tensor product of n copies 
of f). Then the bosonic Fock space over f): 

T = T% = e n >o ® n s f) 

is the space of sequences (p = (v?n)n>o> with (po £ C, (p n £ <8>™h, and with an inner product 
defined by 

n>0 

where (ip n , ifi n ) denotes the inner product of <8>™f). The vector $7 = (1, 0, . . .) G T is called the 
vacuum. By C T we denote the dense subspace of vectors ip for which ip n = 0, for all but 
finitely many n. The number operator Nf in T is defined by (Nf(p) n = nip n . 
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B.l Creation- and Annihilation Operators 

The creation operator a*(h), h £ f), on C T is defined by (a*(h)<p)o = and 

(a*(h)ip) n = \/nS n (h ® <p n -i) 

where S n £ L(<g> n f)) denotes the orthogonal projection onto the symmetric subspace (8>™f) C 
(g> n f). The annihilation operator a(h) is the adjoint of a*(h) restricted to !F^ n . Creation- and 
annihilation operators satisfy the canonical commutation relations (CCR) 

[a(g),a*(h)] = (g,h), [a*(g), a*(h)} = 0. 

In particular [a(h),a*(h)] = \\h\\ 2 . From the definition of a*(h) it is easy to see that 

\\a#(h)(N + l)~ 1 / 2 \\<\\h\\. (20) 

In the case where f) is the one-photon Hilbert space, L 2 (K 3 ;C 2 ), the annihilation and 
creation operators can be expressed in terms of the operator- valued distributions a\(k) and 
a* x (k) by 

a *( h ) = E / hx(k)a* x (k)d 3 k. 

Setting G Xi \(k) = \k\~ 1 / 2 E\(k)x{\k\<A} e ~ tk x i the quantized vector potential A{x) can be de- 
fined as A(x) = a(G x ) + a*(G x ). 



B.2 Second Quantization 

Suppose b is a bounded operator on f) and ||6|| < 1. The operator F(b) : — > -^(f)) is 
defined by 

r(6)n = n 

r(6)t<g£f) = 6®...® 6. 

Clearly ||r(6)|| < 1. From the definition of a*(/i) it easily follows that 

T(b)a*(h) = a*{bh)T{b) (21) 

T{b)a(b*h) = a{h)T(b), (22) 

and hence that T(b)a(h) = a(bh)T(b) if b*b = 1. 
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If b : D(b) C H — > TL is self-adjoint, then dr(6) in .F(h) is defined by 

dr(6)0 = 

n 

&Y{b)\ ® n s D(b) = 1^__1) 

i =1 j-l n-j 

and by linear extension. dr(6) is essentially self-adjoint and, denoting the closure by dr(6) 
as well, T(e lb ) = e* dr W . One example is the number operator Nf = dr(l), another one, for 
h = L 2 (M 3 ;C 2 ), is the field energy 

H f = dT(\k\)= [ \k\al(k)a x (k)d 3 k. 

A=l,2 ^ 

B.3 The Identification Operator 1 : T ® T ^ T 

In the proof of Theorem|f)]an important role is played by the identification operator / : T®T — > 
T defined by 

I(<^(g>Q) = ip 

lip® a* {hi) ■ ■ ■ a*(h n )£l = a*(hi) ■ ■ ■ a*(h n )<p, ip G J^fin, 

and linear extension to ^"fin^-^fin- This operator is unbounded. We often use the commutation 
relation 

a(h)I = I(a{h) <g> 1 + 1 <g> a(h)), (23) 
which is in contrast to a*(h)I = I(a*(h) 1) = 1(1 (g) a*(h)). 
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